Abstract 

In terms of linearized Gross-Pitaevskii equation we have studied the 
process of sound emission arises from a supersonic particle motion in a 
Bose-condensed gas. By analogy with the method used for description of 
Vavilov-Cherenkov phenomenon, we have found a friction work created by 
the particle generated condensate polarization. For comparison we have 
found radiation intensity of excitations. Both methods gives the same 
result. 
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1 Introduction 

The recent experiments |6j of Bose-Einstein condensation (BEC) in trapped 
alkali- metal atomic gases have stimulated a great interest on the problem 
One of the reasons of such an interest is the fact that nonlinear Shrodinger equa- 
tion (Gross-Pitaevskii equation) is well suited for describing the fundamental 
characteristics of an interacting Bose-gas at zero temperature and also appears 
in other physical problems: theory of superconductivity, theory of elementary 
particles, nonlinear optics [|l0| etc. 

Examination of fast particle's energy loss is one of the popular experimental 
technique in condensed matter physics. Such experiments also can give an 
additional information about properties of an ultracold Bose-gas. 

It is known, when a particle moves in a matter with a velocity v above the 
supersonic speed c which equals to the speed of excitations, the matter begins 
to radiate. In electrodynamics this phenomenon names Vavilov- Cherenkov ra- 
diation , in aerodynamics such an effect appears from a supersonic motion of 
a body, in condensed matter physics — this is a polaron problem Q]. 

In the paper we use linearized Gross-Pitaevsky equation to study the process 
of radiation of the excitations in a Bose-condensate, arises from a supersonic 
particle motion. Under this consideration, the Bose-condensate excitations are 
described as the excitations of a classic nonrelativistic scalar field, same as radi- 
ation of a classic electromagnetic field in a continuous matter. As it was shown 
[||, the GPE can be reduced to the system of hydrodynamic equations. But 
this approach adequately describes only the long-wave part of the spectrum, 
which corresponds to the low energy levels of elementary excitations <C /x. 
As it will be shown below, the maximum of the radiation is in the region 
Wmax = 2/i(ti/c)-y/(w/c)^ — 1. So that, direct investigation of the GPE gives 
more detailed description of the problem. Note that to find an energy loss of a 
particle it is enough to solve the problem by the methods of quantum mechanics 
using "Fermi's golden rule." The result will be the same. 
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2 Particle driven condensate polarization 



The state of the condensate (with a fixed chemical potential /i), is described 
by the macroscopic wave-function 5* (r, t) , which satisfies the GPE equation 
(fi=l), 



. d 



■^{v,t) = (^|- + ^^(r,<)-^^ *(r,i)+ j driC/(r-ri)**(ri,t)^(ri,t)^(r,i), 



The interaction energy of Bose-gas atoms [/ (r — ri) = Uq5 (r — ri) and poten- 
tial energy of its interaction with the straightly and uniformly moving particle 
V {r.t) = Vbe**(5(r — vt) ,s = +0. In the last expression we use an adiabatic 
switching of interaction to exclude automatically advanced solutions . Wave- 
function of the ground state of a static gas without a particle corresponding 
to zero energy (counted from the chemical potential) doesn't depend on the 
coordinates and equals ^'o — fi/Uo. In the presence of a particle condensate 
becomes heterogeneous — "polarizes." The "polarization" , in general, depends 
on time ip (r,t) — (r,t) — ^o, iu linear approximation it satisfies inhomoge- 
neous equation: 

z^V (r, = (l^ + m) (r, t) + fir (r, t) + be^'S (r - vt) , (1) 
6 = ^0*0, [5] = p/Vi],s = +0. 

Let's solve the equation (|l|) using Fourier's method. We present the potential 
and the 5— function as a plane waves decomposition, 



V'(r,t)= / ^^e^''-— *+^V(k), (2) 



6{r- vt) e"* = / 7^6*'^'"-*"*+"*, lo = kv. (3) 



(2^ 

Substituting (||) and (^) into equation (|l|), we get 

{u; + ^s-Oi^i^)^^^r{-K) + b, ^ = + (4) 

2m 

Combining this equation with complex conjugated one, we find 



^(k)=6 ^^^^^ (5) 



2m \ 2m / m 

The poles of (||) describes the Bogolyubov's oscillations spectrum e (fc) of the 
Bose-condensate. Substitution of (|^) into (|^) gives an important expression for 
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the field of condensate "polarization" , 

With a slow particle motion the field decays exponentially at a large distance 
from the particle. The poles of (j^) can lead to radiation, 

kv = ±e (fc) 

Thus, for appearance of radiation in a sound spectrum region (e (fc) — ck) we 
have the condition, 

fct;cos0o = cfc cos^o = ~- 

V 

It gives the equation of the front of sound radiation. With a particle motion 
parallel to z axis, we have 

z — vt — asin6o, r = acos9o, < a < oo. (8) 

If we place a center of spherical coordinates at the point of the particle then 
there is no field at right of the front in the region < 6 < 6q + ^tt. Properly 
speaking if we consider the dispersion of sound speed then high-frequency field 
in the region is not zero, but only the low-frequency part of the spectrum 

C (fc) = ^-7-^ = \/ + -^fc2 < y 

k V 4m^ 

gives contribution to radiation. The dispersion of the sound phase velocity 
c (fc) plays the same role as the dispersion of inductivity, cutting the radiation 
frequency on top, 

C(fcmax) = w, (9) 
t^max = fcmaxW = £ (fcmax) = 2mV\/ - (10) 



3 Fast particle energy loss 

In first order of perturbation theory, the potential energy of the particle in the 
field of condensate "polarization" is 

6E = j dv-^* (r, t) y (r - R) (r, t) , 

= 6*o + 2&ReV'(R), 

where R(t) — is the radius vector of the moving particle. The energy produces 
the drag force which acts on the particle from its created field 



dR 



-5E 



-26 Re 



rfR 
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This results in particle energy loss 



E = Fv = -26vRe 



that retires to infinity as a form of condensate excitations. Using the expression 
we find the radiation intensity 



-E = 2b^ Re 



(ikv) ■ 



2m 



To integrate the expression let's use Sokhotsky formula 

^ =P^^-^[J(c.-e) + 5(c. + e)] 



(w + isY - 

The principal value of the integral gives no contribution into energy loss. Thus, 
/= k^dk J dcose (^■^+^^ [S {lu - e) + 5 {uj + e)]. (11) 

-1 

Transforming ([ll| ) using a new variable oj = kv cos 9, we get the very simple 
formula, 



/ = 



k^dk. 



(12) 



With the help of the formula which is inversed to (j^) , 

k^ = 2m(-/i + ^J'iJF+7^), 
we have an expression for the intensity in form of a spectrum decomposition 



2ttv 



s\ 1 



de. 



(13) 



4 Radiation intensity 

Now, it is interesting to directly evaluate the radiation intensity, integrating the 
energy flow Q by all directions |1| , 









Re 




4^ 
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. d 

where P = — iV — is the operator of momentum, i— — the operator of energy. 

Whole energy flow through the surface of the cylinder with radius i? — > cx) 
around the particle trajectory is 

I = [ QdS = 27:RRe [ dz ( ^-^A t^^,p (14) 



To evaluate the integral let's find form of the field (Q) in the wave zone in 
cylindrical coordinates 



{lo + isY — £^ (fc) ' 



(15) 



Introducing the "retarded" time t' = t — —, and expressing the polar angle 

V 

integral (by variable ip) using the definition of Bessel function 

r2TT 

gZxcosv.^^^27rJo(x), 



we find 

oo oo 

^(f,i) = / k^dk±Jo{k±R) / (iu;e-*("+*")* - 

(27r) Am^{uj + is)^-4mfi{kl + ^)-{kl + 

to get more compact expression let's factorize the denominator 



where 



(27r) w 



duje 



F, F ^ J 2k^dk^Ja ik±R) 





{kl - kl) {kl + kl) - isLO 

(16) 



kl = -{ 2fim + ^ ) + 2mv/tj2 + ^2 



kl = 2fim + ^ + 2my/LU^+Jj? 



With small frequencies the positive value of k\ ~ uP' have a maximum at 
u! = \ {mv^)'^ — /i^ and becomes zero at Wmax = 2mv\/ v"^ — c. The integral 



(w) can be represented in a following form 



r,2 , 



i^(t^) = / dk^Jo(k^R)- 



2muj 



2fci 



- (yfei + iujsy K + 



(17) 
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The last term in the square brackets has imaginary poles, it gives the expo- 
nentially small contribution to the integral at large distances from the particle 
trajectory {R — > oo) and can be neglected. The expression ( pT[ ) can be trans- 
formed to integral by whole real axis Q . 

Let's examine the integral along the closed contour 

J=i| f{z')H^'^iz)zdz (18) 
c 

where Hq^"^ (z-R) — is Hankel function, analytical in an upper half-plane of 
complex variable z which have the following asymptotic behavior on the infinity 

H^q\z)k X —e''-''''\ (19) 

V TTZ 

The path C consists of the real axis and the infinite semicircle in an upper- 
halfplane. The integral along the last contour is exponentially small. So we can 
write 



oo 



\ J dzH^^^ (zR) zf (z2) , J+ = IJ dzH^^^ (zR) zf (z') 

-oo 

(20) 

Let's transform first integral to a new variable z = pe*''' 



J-=-\j dp4'\e^-pR)pf{p'). (21) 



Taking into account ( pO| ) we have for sum ( pl| ) 

oc oo 

J=IJ dppf {p') [/Je^''^--/4 + yAe-P«+-/4] ^ J dppf (p2) j„ i^pR) . 



On the other hand the integral on the closed contour (^8|) equals the sum of 
residues of the poles Zn in upper-halfplane 

J = 2^iY, Itj'^) {znR) z„f (zl) (22) 

res 

Applying (|2^) to integral (|l^) and remembering that integrand have one pole 
at < Wmax close to the real axis in the upper-halfplane equal to kj_ = 
(fci signw -I- is) we have 
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F {u;) = 2mH^'^ {p) \ + "A^^"" = vrzi/^^^ {p) D (a;) , (23) 
D{lo)^^ f =1+ ^ =, fciE signer + (24) 



Substituting ( p3[ ) into ( |16| ) and taking into account that at |a;| > Wmax parameter 
fci becomes imaginary and Hankel function is exponentially small, we find ■0 (r, t) 
in the wave zone as a form of single integral by frequency 

V'(r,t) = -^ j d^p-(-+-)*'ijW(p)i^(c.). (25) 
Let's use this expression to find (|lj). Using (^5|), ([l9|), (24) we have 

_^4;^ = _!;!!^ / duj(^siqnoj\e-^^-+^^^''H^^\p)D{u:) (26) 



and 



= -i^ / d.(u.)e-(--)*'<)(p)i.(u.) (27) 



Substituting this expressions into (|lj) and integrating by z with a formula 

-OO 

and using expHcit form (|l9|) we have 



\8ttv J J m t:\p\ 



/ du}\u!\D (uj)^ ~ / ujduj ( 1 , ^ 

8ttv J ^ ^ ^ ' 2tiv J I JJ?T^ 



The last expression is identically coincides with (^^ . We see that radiation in- 
tensity is monotonically grows with the frequency. To evaluate integral radiation 
intensity it is the simplest to integrate expression (IlJ 
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So, we have the exact agreement of the expressions for the radiation spectrum 
of the bose-condensate excitations generated by the fast particle received by 
the drag method and the method of the radiated energy flow evaluation. It is 
possible to show the "Fermi's golden rule" gives the same result. Thus near the 
threshold the radiation intensity grows as {v"^ — (?Y ■, but at high velocity — it 
is as u^. If in the case of Cherenkov effect the particle velocity is limited from 
the top, then for the nonrelativistic Bose-gas there is no such a limit and with 
f ^ c the radiation intensity can be very big. 

We thank RBRF and INTAS for finance support. Grants INTAS-97-0972 
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